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Abstract. We report the first results of a study of a sample of 20 galaxy mergers/interacting systems, using the VIMOS and 
PMAS integral field spectrographs. For each object, we extracted the integrated spectrum of the central regions and analyzed 
the ionization state using classical diagnostic diagrams l*Veilleux"& Osterbro ckll987l) . There is evidence of AGN ionization in 
4 of the objects, i.e. 20% of the sample, a considerably higher fraction than found in previous studies (~4%). An increase of the 
nuclear activity with the far infrared temperature may account for the difference. 



Key words. 



1. Introduction 

Interactions and mergers can strongly affect galaxy e volu- 
tion. Since the early stud ies bv lToomre & Toomre I (11972) and 
Larson & Tinslev ( 1978), it is generally believed that interac- 
tions can trigger starburst processes in galaxies, via tidal ef- 
fects. The fact that interacting galaxi es are more lu minous in 
the far in frared (FIR) (e. g. ISanders & Mir abel 1996) and the 
radio (e.g.lHummel 198 Tb, and that they have larger Ha equiv- 
alent widths (e.g. lKennictitt et alJl987l) supports this hypothe- 
sis. The star formation is found to be stronger in the nucleus, 
although it is also enhanced in the disk. A widely accepted 
idea, supported by simulations, is that galaxy interactions in- 
duce flows of gas from the outer parts of the galaxy and/or 
the companion into the inner regions throu gh loss of the angu- 
lar m omentum induced by tidal forces (e.g. lMihos & Hern auist 
1996). 

It has long been suggested that strong interactions and 
mergers may trigger nuclear activity. Many authors have 
found that quasar host galaxies show di storted morpholo- 
gies, rem i niscent of past merging events (e.g.lMcLeod & Rieke 
19971 ISanchez & Gonzalez-Serrano 
2003; Jahnke et al.l 120041). However, other authors found no 



evidence of past interactions (e.g. lDunlopetal]l2003l) . It is 
not clear if the fraction of AGN host galaxies undergoing in- 
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teractions is really larger than the fraction of inactive galax- 
ies undergoing interactions, or how the timescales of the 
morphological relaxation and the AGN ignition relate (e.g. 
Sanch ez et al . 2004). However, it is clear that some AGN, if 
not all, are gene rated b y mergers/i nteractio ns betwe en galaxies 
JCanalizo & StocktonfeOOlUCattaneo & et all2004l) . 

The fraction of AGNs in merging systems increases with 
the FIR luminosity, being ~30% for the e xtreme case of 
the ul traluminous infrared galaxies (ULIRGs) (Bushouse et al. 
2002k with FIR luminosity L F!R > 10 12 Lm. ISanders etalJ 
( 1988) suggested an evolutionary scenario, in which major 
mergers may allow the infall of gas towards the nuclear re- 
gions, triggering nuclear star formation and igniting an AGN. 
Both the AGN and the star formation heat th e dust , and the ob- 
ject is observed as a ULIRG. IColina e t al. (2001) proposed a 
modification of this scenario in which cool ULIRGs would be 
generated during the merging of two or more low-mass sub- 
U spirals, and warm ULIRGs in the merging of an L* galaxy 
with a substantially less massive galaxy or during the merging 
of two intermediate-mass spirals. Under this hypothesis only 
warm ULIRGs would evolve/host an AGN, while cool ULIRGs 
would be heated by star formation. 

At lower luminosities, L FIR ~ 1O 1O_I1 L , there 
is no evidence of an excess of AGNs among inter- 
acting/mer ging systems compared to i solated galax- 
ies ( e. g. iDaharil 1 19851 iKeelet alJ 1 19851 iBushouse et all 
1988; Sekiguchi & Wolstencro ftl 1 19921 Liu & K ennicutt 
1995 



Sekiguchi & Wolstencroft 

IPastoriza et alJ 1 19991; iDonzelli & Pastorizal bOOCt 
Ber gvall et al.l 120031) . This deficiency could be partially 
explained by the presence of strong circumnuclear star 
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formation that can blur the AGN signatures JPastoriza et alJ 
ll999tlPonzelli & Pastorizal EoOO). There is significantly more 
nuclear star format ion in interacting ga laxies, compared to 
isolated ones (e.g. iBergvall et al.l 120031) . T his nuclear star 
forma t ion increases with the FIR temperature ( Bush ouse et alJ 
1988; Berevall et al. 2003). If the scenario proposed by 
Coli na et al.l (1200 11) is valid for lower-luminosity objects, the 
nuclear SFR and the probability of generating an AGN c ould 
be larg er in warm FIR sources than in coo l ones. lde Griip et al~l 
i 19851) and Ivan den Broek et al~l J 1 99 lit already noticed an 
increase of the number of AGNs with the FIR temperature, 
using slit spectroscopy. 

Most previous studies of this topic were based on samples 
of in teracting/mer ging galaxies selected by mo rphology (e.g. 
iDonzelli & Pastorizall997UBergvall et alJl2003l) . without con- 
sidering the FIR properties. As a consequence, these samples 

contain mainly cool FIR sources (log I( )(/60/j/ /lOO/a) 0.4). We 

selected a sample of 20 IRA S sources with evidence of interac- 
tions/major merging ( S anchez et alJ2004bl , paper I), with sim- 
ilar FIR luminosities to previous samples, Lfir ~ 10 10_11 Lo, 
but warmer (logio(/60/i//loop) ~ -0.25). Since the use of 
slit spectroscopy can bias the results (paper I), we obtained 
integral-field spectra of all them to explore the connection be- 
tween galaxy interactions and nuclear activity. 

2. Observations and data reduction 

The spectra were taken with the integral-field units (IFUs) 
of VIMOS JLeFevre et all l2003fc . 17 objects, and PMAS 
jRoth et al.ll2000l) . 3 objects. The VIMOS observations were 
made on the nights of the 9th, 10th, 22th, 23th, 24th and 26th 
of April 2003 at the VLT/UT3, in clear sky conditions, with 
an average seeing -1.6". We obtained, for each object, 3 ex- 
posures of 300s using each of the low-resolution red (LRr) 
and blue (LRb) grisms. With this setup we cover the wave- 
length ranges 6800-9500 A (LRr) and 3900-6800 A (LRb) with 
resolutions of 7.14 A/pixel (LRr) and 5.35 A/pixel (LRb), re- 
spectively. The VIMOS IFU contains a focal-plane array of 
6400 lenslets each of size 0.67"x0.67" (in the low-resolution 
mode) distributed on a regular grid. The final field-of-view was 
54"x54", which allows us to cover the different components 
within the IRAS aperture (~1')- In two cases (IRAS 10219- 
2828 and IRAS 12110-3412) we obtained data for only one 
grism (LRr and LRb, respectively). 

The PMAS data were taken during the nights of the 28th 
of April, and 5th and 8th of August 2003 at the 3.5-m tele- 
scope of the Calar Alto Observatory. The observations were 
made in poor weather conditions as a back-up project for these 
nights. We used the V300 grating which covers the wavelength 
range ~4000 - 7700 A with a resolution of 3.3 A/pix. The exact 
wavelength range depends on the rotation angle selected for the 
grating, which changes from object to object. The PMAS focal 
plane contains an array of 256 lenslets each of size 0.5"x0.5" 
(in the standard configuration) distributed on a square grid. The 
field-of-view is 8"x8", i.e. too small to cover the IRAS aper- 
ture with a single pointing. Different mosaic patterns were used 
for each object to cover the central regions and different regions 
of interest (e.g., tidal tails). Due to the variable atmospheric 
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Fig. 1. Intensity map of the continnum emission of IRAS 
13477-4848 at wavelengths -6477-6591 A, extracted from 
the VIMOS datacube. The contours show the emission in 
Ha+fNII] on a logarithm scale, obtained by subtracting the 
continuum emission from that at wavelengths centred on the 
emission lines (~6612-6670 A). 



conditions different exposure times were used for each object 
and each night, to obtain similar signal-to-noise levels in all the 
frames. The exposure times ranged from 1200 s per pointing for 
IRAS 16365+4202 (05/08/03) to 5400 s for IRAS 12110+1624 
(28/04/03). 

Both the VIMOS and PMAS data were reduced using 
the P3d reduction package, written in IDL llBeckerl 12002). 
Although initially written to reduce PMAS data, it was slightly 
modified to reduce da ta from other IFUs, such as VIMOS (e.g. 
ISanchez et alJl2004ah . The expected locations of the spectra 
were traced on a continuum-lamp exposure obtained before 
each target exposure. After bias subtraction, the target spectra 
were extracted by adding the signal from the 5 pixels around 
the central traced pixel. The continuum-lamp exposure was 
also used to determine the fiber-to-fiber response of the instru- 
ment at each particular wavelength (the so-called flat-spectra). 
These flat-spectra were used to homogenize the response of all 
the fibers. Wavelength calibration was carried out via observa- 
tion of emission-line lamps, also obtained before the science 
exposures. The accuracy of the wavelength calibration was 
checked using sky emission lines, and has standard deviation 
1.2 A (VIMOS-LRr), 1.0 A (VIMOS-LRb) and 0.5 A (PMAS). 
Different exposures taken at the same pointing were combined 
using IRAF tasks. A 3cr clipping algorithm was used to remove 
possible cosmic rays. The different PMAS pointings were mo- 
saiced into a singl e frame, using the Euro3D data format, as 
described in (Sanchez et al. 2004a). 



S.F.Sanchez et al.: The merging/AGN connection 



3 



Further analysis was carried out using E3D ( ISanchezl2004 . 
and our own fitting routines (Sanchez et al., in preparation). 
An average sky spectrum, taken from areas uncontaminated by 
target spectra, was subtracted from the target spectra. 

3. Analysis and results 

3.1. Extraction of the nuclear spectra 

Figure [2 shows an example of the reduced data: a continuum- 
intensity map of IRAS 13477-4848 at wavelengths -6477- 
6591 A, extracted from the VIMOS datacube. For each ob- 
ject we extracted the spectrum of the central region, by co- 
adding the spectra within a box aperture of size 2"x2" (ap- 
proximately the central kiloparsec, at the redshift of these ob- 
jects). Whenever the intensity peak of the continuum emission 
was not coincident with that of the emission line, we centered 
our aperture on the latter to increase the signal-to-noise ratio. In 
most cases both peaks were coincident within a couple of arc- 
sec. This is illustrated in Fig.^ where the Ha+[NII] intensity 
contours are superimposed on a continuum image. 

Six of the 20 objects have two well-separated components, 
with two clear intensity peaks. The spectra of four of them 
show emission lines in both components. Emission lines are 
absent in the remaining two. The extracted spectra of the nu- 
clear regions of the two components have similar characteris- 
tics. Hereafter we restrict our discussion to the properties of the 
spectrum of the brighter component. 

3.2. Ionization Mechanism 

The final dataset consists of the spectra of the nuclear re- 
gions of 20 interacting/merging galaxies. Of these, 18 cover 
a wavelength range which includes some of the most important 
emission lines for ionization diagnostics: i.e., H/3, [OIII]/l4959, 
[OIII]i5007, Ha, [NII]^6584, , [NII]/16548, [SIIJA6717 and 
[SII]/16731. A variety of other emission lines were detected 
in many of these spectra. A detailed description of the spec- 
trum of each object will be presented elsewhere (Sanchez et 
al., in preparation). We fitted the line profiles to derive the in- 
tegrated flux of each emission line, using our own software. 
A single gaussian was fitted to each line, using a polynomial 
function to characterize the continuum. The Ha- [Nil] system 
and the [SII] doublet were blended in the LRr VIMOS data. 
For proper deblending, the minimun line-width was fixed to 
be the width measured in isolated sky emission lines (i.e., the 
instrumental width). Our software allows definition of emis- 
sion line systems, i.e., a kinematically coupled set of emission 
lines with the same width and heliocentric velocity. This was 
essential for accurate deblending of the lines. After deriving 
the flux of each line, the [01111^5007/^, [NII]^6583/Ha and 
[SII]/167 17+673 1/Ha line ratios were calculated. 

Figure [5] shows the basic diagnostic diagram, 
[OHIU 5007/Hg vs. [NII]/l6583/Ha JVeilleux & Osterbrockl 
1987), for the 18 objects observed with spectral informa- 
tion in both wavelength regions (solid squares). The open 
squares near the edge of the plot represent two objects (IRAS 
10219-2828 and IRAS 12110-3412) with data only for one 



VIMOS grism (LR-red or LR-blue). The solid lines show the 
classical separation between ionization due to star formation 
(bottom-left), AGNs (top-right) and LINERs (bottom-right). 
The boundary curve between AGNs and star forming regions 
is taken from IVeilleux & Osterbrockl Jl987t) . A separation 
between AGNs and LINERs at a ratio of [OIII]/H/3 = 3 
was assumed, following earlier results in the literature (e.g. 
iDonzelli & PastorizafeOOollGarcfa-Lorenzo et all200ll) . 

The majority of the objects in our sample (14 out of 20) 
have nuclear emission features characteristic of star forming 
regions. However, for a significant fraction of the objects (4 
out of 20, 20%), the nuclear emission is dominated by an 
AGN: IRAS 16330-6820, IRAS 12193-3942, IRAS 13031- 
5717 and IRAS 16365+4202. We include the lcr error bars 
for these objects in Fig. 2. As a check, we repeated the anal- 
ysis using the [OIII]/l5007/H/3 vs. [SII]^6717+6731/Ha di- 
agnostic diagram. We find the same distribution, reinforc- 
ing the result from Fig 2. As expected from ionization mod- 
els re.g. lOsterbrockl ll989). there is a correlation between the 
[NIIR6583/Haand the [SII]^6717+6731/Haline ratios in our 
sample, with a dispersion dominated by measurement errors. 
The errors in the [NII]/l6583/Ha ratio are dominated by resid- 
ual errors in the line-deblending procedure, while those of 
[SII]/l6717+6731/Ha (less affected by deblending) are dom- 
inated by the lower signal-to-noise ratio of [SII]/16717+6731. 

For the two objects with only one diagnostic line ratio mea- 
sured, the nature of the ionization is uncertain , although one of 
them, IRAS 121 10-3412, is clearly not an AGN. 

4. Discussion and Conclusions 

As noted above, most studies of the ionization mechanisms in 
the nuclear regions of merging/interacting galaxies are based 
on studies of cool FIR sources. These studies found a small 
fract ion of objects with AGN-like ionization, typically ~4% 
fe.g. IPastoriza et aljfl99^ IBergvall et alJl2003h . By contrast, 
we find a significantly higher fraction of AGNs (~20%), in a 
sample of warmer FIR merging/interacting galaxies. This may 
indicate that the scenario proposed bv lColina et alJ J200ll) (that 
only warm ones host an AGN), is valid also at lower IR lumi- 
nosities LflR. 

Our current sample corresponds to only ~6% of all interact- 
ing syste ms with measured r edshift in the IRAS point-source 
sample JSanchez et alj r2004b). How general is our result? To 
addr ess this question we compared our sample with that of 
IBer gvall et ail i200A As already mentioned, the two sam- 
ples cover the same range of FIR luminosities (indistinguish- 
able, according to a Kolmogorov-Smirnov test). However, the 
ranges of FIR temperatures are different, with our sample being 
warmer (the null hyphothesis, that both samples were drawn 
from the same parent population, can be rejected at the 0.13 % 
level). If we add the objects from Bergvall et al. (2003) within 
the IRAS point source catalogue Jjointlras Sciencdl 19941) to 
our sample, we end up with 77 objects. Splitting this sample 
into cool (feon/fiooii < -0.3) and warm (/60jii//loo/j > —0.3) 
FIR sources, we found that only 4% (2 out of 50) of the cool 
objects harbour an AGN, in compared with 15% (4 out of 27) 
of the warm ones. That is, we still find that the fraction of AGN 
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Fig. 2. [OIII],15007/H V 3 line ratio as a function of 
[NII]/K)583/HaT for the 18 objects with both ratios mea- 
sured (solid squares). The open squares indicate the two 
objects with data only for one VIMOS grism (LR-red or 
LR-blue). They have been plotted with an arbitrary value 
of the un-measured line ratio (i.e. at the edge of the figure). 
The solid lines show the separation between ionization due 
to star formation (bottom-left), AGN (top-right) and LINERs 
(bottom-rig ht). The approximate locus of star-forming objects 
reported bv lVeilleux & Osterbrockl dl987l) . is indicated with a 
dashed line. 

increases with the FIR temperature. We found similar results 
when com paring with other samples (e.g. lDonzelli & Pastorizal 
119971120001) . 

The probability of igniting an AGN seems to be higher 
in warm FIR sources, over a larg e range of FIR luminosities. 
According to IColina et alJ (|2001), these objects are the result 
of a merger either between a large galaxy and a small galaxy 
or between two galaxies of similar size. A possible explanation 
of this result could be that apart from the fuel supply provided 
by the interaction/merging, the presence of a massive black- 
hole is required to ignite an AGN. Massive black-holes are 
found in mass ive galaxies, due to th e black-hole/bulge-mass 
relatio n (e.g., iMagorrian et alJ fl~998t iKormendv & Gebhardt 
2001 ). Therefore, only galaxy mergers with one progenitor har- 
bouring a massive black hole (i.e. a massive galaxy) may (re- 
ignite an AGN, and those are the mergers that generate a warm 
FIR source. 
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